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ABSTRACT 

Fatigue  reliability  calculations  during  iterations  of  the  design  process  have  been  demonstrated  as  practical  using  new 
analytical  and  numerical  methods  based  on  analytical  cumulative-damage  reliability  tools  (a-functions).  Effects  of 
fatigue  strength  and  usage  distribution  may  be  studied  in  the  presence  of  load  variation.  However,  to  date,  the 
development  of  the  a-functions  has  been  centered  on  use  of  the  Palmgren-Miner  cycle-ratio  summation  rule  for 
cumulative  damage,  also  known  as  the  Linear  Damage  Rule  (LDR).  This  paper  extends  the  a-function  toolbox  for 
use  in  determining  reliability  based  on  the  Manson-Halford  Double  Linear  Damage  Rule  (DLDR).  Examples 
demonstrate  use  of  the  method  and  charts  are  provided  to  illustrate  the  sensitivity  of  the  cumulative-damage 
reliability  problem  to  load  and  strength  variations.  For  a  given  reliability,  probabilistic  DLDR  methods  demonstrate 
significant  reductions  in  life  when  compared  to  LDR  methods.  The  importance  of  understanding  and  applying  the 
appropriate  LDR,  DLDR,  or  other  material  characterization  is  apparent. 


NOTATION 

A  material  constant  (log-log,  intercept) 

a(»,«)  Benton  a-function  defined  in  [5] 

a*  (•,«,•)  generalized  a-function  defined  in  [7] 
a*DLDm  (•-•/•/•/•)  generalized  DLDR-based  phase  I 
a-function  defined  herein 

a*DLDmi  (•/•/•/•-•)  generalized  DLDR-based  phase  II 
a-function  defined  herein 
B  material  constant  (log-log,  slope) 

C  material  constant  (log-log,  cycles  at 

endurance) 

COVE  endurance  limit  coefficient  of  variation 
COVs  load  amplitude  coefficient  of  variation 

D  accumulated  fatigue  damage 

D,_n  accumulated  fatigue  damage  threshold 

between  DLDR  phases  I  and  II 
Di  DLDR  phase  I  damage  fraction 

Dn  DLDR  phase  II  damage  fraction 

E  endurance  limit 

Ew  working  endurance  limit,  pE(l  —  3  COVE ) 

E§  critical  endurance  limit  for  D  =  S 

exp(*)  exponential  function 

fE  (•)  endurance  limit  probability  density 

function 

/En(«)  E-n  joint  probability  density  function 

fn  (•)  damaging  cycles  probability  density 

function 
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cyE(*)  damage  rate  as  a  function  of  endurance 
limit 

QEn  (***)  accumulated  damage  as  a  function  of  E 
and  n 

t  load  cycle  index 

j  component  index 

k  material  constant  (factor) 

m  material  constant  (exponent) 

m0  DLDR  a-  (and  ^-(function  parameter 

(exponent) 

N  number  of  cycles  to  failure 

Nrej-  DLDR  reference  number  of  cycles 
N*  number  of  cycles  at  DLDR  break  point 

N[  number  of  cycles  within  DLDR  phase  I 

Nn  number  of  cycles  within  DLDR  phase  II 

n  number  of  applied  load  cycles 

nI  number  of  applied  phase  I  load  cycles 

nn  number  of  applied  phase  II  load  cycles 

P  variable  used  for  DLDR  (factor) 

p  load  condition  index 

Q  variable  used  for  DLDR  (exponent) 

R  component  reliability 

S  oscillatory  load  amplitude 

S95  95th  percentile  oscillatory  load  amplitude 

(estimated  “top  of  scatter”  load) 
v  generalized  a-function  parameter 

x  dummy  variable  used  in  Gauss-Laguerre 

quadrature 

y  DLDR  a-  (and  &-)function  parameter 

(factor) 

z  standardized  normal  load  amplitude 

z0  standardized  normal  damaging  load 

amplitude 
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S  accumulated  fatigue  damage  at  failure, 

phase  I  damage  fraction  at  phase  I-II 
threshold,  or  phase  II  damage  fraction  at 
failure  (typically  5=1) 

SA  material  constant  tolerance  band  (log-log, 

intercept) 

/r  mean  (of  load  amplitudes) 

(UE  mean  endurance  limit 

/xs  mean  load  amplitude 

<f  load  condition  usage  spectrum  weighting 

factor 

a  standard  deviation  (of  load  amplitudes) 

</>(•)  standardized  normal  probability  density 

function 

[  ]  vector 

[]r  vector  transpose 


Readers  wishing  to  understand  differences  between 
Palgrem-Miner  LDR  and  Manson-Halford  DLDR 
theories  and  various  other  cumulative  fatigue  theories 
are  referred  to  the  literature  survey  in  [1 1], 

BACKGROUND 

ASTM  Standard  E  739  [12]  equation  (2a)  presents  a 
common  characterization  of  fatigue  strength,  namely: 

log N  -  A  +  B  (logS)  (1) 

which  is  applicable  for  values  of  S  above  the 
endurance  limit.  This  “log-log”  representation  of 
fatigue  strength  is  typically  plotted  with  tolerance 
bands  as: 

log  N  =  A  +  B  (log  5)  ±  8a  (2) 


INTRODUCTION 

Recent  advances  in  helicopter  fatigue  engineering 
methods  have  demonstrated  great  promise  for 
providing  efficient  calculations  to  determine  the 
fatigue  reliability  of  helicopter  dynamic  components 
(see  [1], [2], [3], [4], [5], [6],  and  [7]  for  example). 

Once  mature,  such  methods  will  enable  helicopter 
materiel  developers  and  original  equipment 
manufacturers  to  design  more  robust  equipment 
while  reducing  the  life-cycle  costs  inherent  to  system 
safety  risk  management,  airworthiness  qualification, 
and  continued  airworthiness  processes. 


Let  Nel  =  10c  represent  the  number  of  cycles  at  the 
endurance  limit,  where  6  <  C  <  10  and  the 
logarithm  of  the  corresponding  endurance  limit  is 

l°g  E  =  c-^  (3) 


As  explained  in  [7],  the  log-log  strength  curve  may 
be  rearranged  to  calculate  the  damage  fraction  for  the 
i-th  cycle  applied  to  the  j- th  component  as 


if  Si  >  Ej 
otherwise 


(4) 


In  [5]  and  [7],  the  author  has  developed  and 
demonstrated  new  cumulative-damage  reliability 
tools  (so  called,  “a-  and  b -functions”)  which  may  be 
used  to  study  effects  of  fatigue  strength  and  usage 
distribution  in  the  presence  of  load  variation. 
However,  to  date,  development  of  a-  and  b-functions 
and  generalized  a-  and  b-functions  has  been  centered 
on  use  of  the  Palmgren-Miner  cycle-ratio  summation 
rule  for  cumulative  damage,  also  known  as  the  Linear 
Damage  Rule  (LDR). 

This  paper  extends  generalized  a-functions  to  allow 
determination  of  fatigue  reliabilities  based  on  an 
alternative  cumulative  damage  theory,  namely  the 
Double  Linear  Damage  Rule  (DLDR)  developed  by 
Manson  and  Halford  (see  [8], [9],  and  [10]  for 
details).  Examples  demonstrate  use  of  the  DLDR  a- 
fimctions.  In  addition,  charts  are  provided  to 
illustrate  the  sensitivity  of  the  DLDR  cumulative- 
damage  reliability  problem  to  load  and  strength 
variations. 


where  k  —  10c  and  m  =  —B  are  material  constants 
and  iVjj  is  the  number  of  cycles  to  failure  for  a  given 
load  St  and  endurance  limit  Ej.  Considering  the 
material  constants  k  and  m  (and  therefore  B  and  C) 
as  invariant  within  the  population  of  components,  the 
logarithm  of  the  endurance  limit  for  the  j-th 
component  is 

log  Ej  =  (5) 

As  derived  in  [7],  the  following  expression  represents 
the  Palmgren-Miner  LDR  damage  accumulated  in 
n  cycles  for  normal-distributed  loads: 


where  fi  is  the  mean  load,  a  is  the  load  standard 
deviation,  and  the  “generalized  a-function”  is  defined 
for  v  <  z0  as 

00 

(z  —  v)m  4>(z)dz  (7) 


a 


with 


(8) 

which  is  the  standardized  normal  probability  density 
function. 

Manson  and  Halford  have  proposed  a  Double  Linear 
Damage  Rule  (DLDR)  composed  of  two  phases  with 
distinct  linear  fatigue  damage  accumulation  patterns. 
Although  the  concept  of  two  phases  originated  from 
an  analogy  to  crack  initiation  and  crack  propagation, 
these  two  phases  are  no  longer  identified  with  a 
definable  physical  process  (see  [10]).  For  the  i-th 
load  level  and  the  j- th  component,  one  may  describe 
the  number  of  cycles  in  each  phase  as  NI  t  j  and  Nu  ij. 

Per  [10],  the  fatigue  damage  accumulated  at  the  end 
of  phase  1  is 

(Nref\02S 

£/-//  =  0.35  (^J  (9) 
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Figure  1:  Double  Linear  Damage  Rule 
Determination  of  DLDR  phases  using  the  N* 
damage-curve  break  point  (adapted  from  [10]). 


where  Nref  is  the  Palmgren-Miner  LDR  number  of 
cycles  for  a  reference  load  level  and  N*  is  the 
particular  life  level  used  to  define  phase  I  with  a 
damage-curve  break  point  occurring  at  damage  D;w/, 
as  shown  in  Figure  1. 

Manson  and  Halford  recommend  that  the  value  of 
Nref  be  selected  lower  than  all  life  levels  in  the 
spectrum.  As  such,  it  is  typically  acceptable  to  select 
Nref  —  103  as  a  lower  bound  for  high-cycle  fatigue 
loads. 


for  the  i-th  load  level  and  the  y -th  component,  the 
number  of  cycles  in  phases  I  and  II  are  calculated  as 

NIU  =  N[j  exp(P/Vt/)  (10) 

and 

N„,i,j  =  Nij  -  Nim  (11) 

respectively,  where  intermediate  variables  Q  and  P 
are  defined  as 


Manson  and  Halford  recommend  in  [8]  that  N*  be 
selected  as  the  highest  life  level  in  the  spectrum  and 
provide  an  example  of  minimal  effects  resulting  from 
alternate  N*  selections.  In  [9]  and  [10],  Manson  and 
Halford  propose  various  more  complex  spectrum- 
dependent  iterative  methods  of  selecting  values  for 
N*.  Notwithstanding  such  recommendations,  the 
probabilistic  nature  of  the  loads  and  spectrum 
inherent  to  present  work  prompts  the  author  to  delay 
consideration  of  D,_;/  spectrum  dependency  by 
pursuing  initial  development  of  DLDR-based 
methods  based  on  selecting  N*  =  NEL  —  10c  (i.e., 
the  highest  damaging  life  level  possible  in  any 
spectrum).  Once  the  author’s  DLDR  methods  are 
established,  readers  may  wish  to  use  the  methods 
derived  herein  to  assess  the  potential  benefits  of  re¬ 
introducing  £>/_//  spectrum  dependency  into  more 
complex  methods. 

Adapting  Manson-Halford  equations  in  [8]  to  avoid 
notation  conflicts  and  facilitate  calculating  damage 
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In 
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InCD,-,,) 


.  ( .  0.65  \ 
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In  £>/-;; 

NrefQ 


(12) 


(13) 


and  where  is  the  LDR  number  of  cycles  to  failure 
for  the  i-th  load  level  and  the  y'-th  component. 

Dashed  lines  in  Figure  1  show  the  results  of  adjusting 
damage  curve  break  points  for  life  levels  other  than 
Nref  and  N*  to  match  the  N,  t  j  =  Ntj  exp(PiVj  yQ) 
equation  based  the  N*  break  point.  Thus,  the  DLDR 
formulation  allows  fatigue  damage  to  be  considered 
linear  in  each  phase  while  closely  approximating 
damage  curves  generated  from  two-load-level  fatigue 
test  results  (see  [8],  [9],  and  [10]  for  these  test  results 
and  other  details  related  to  the  DLDR  formulation). 


- ASTM  Curve  Shape,  A  =  6,  B  =  -2 

. DLDR  Phase  I,  Nref  =  10A3  N*  =  10A6 

- DLDR  Phase  II,  Nref  =  10A3  N*  =  10A6 
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for  LDR.  As  such,  a  new  class  of  generalized 
DLDR-based  phase  I  a-functions  is  defined  herein 
for  v  <  z0  as 

aDLDRl(Z0’m’V’y’mo)  ~  /ic\ 

p  (j^r  (15> 

JZo  exp[y(z  -  u)m°]  * 

where  0(z)  is  the  standardized  normal  probability 
density  function. 

Generalized  DLDR-based  phase  I  a-functions  may  be 
applied  to  phase  I  damage  calculations,  as  follows 

D,j  =  n,  gE}I(Ej)  (16) 

where  Dt  j  is  the  Manson-Halford  DLDR  phase  I 
damage  fraction  and  the  corresponding  phase  I 
damage  rate  gE  i{Ej)  is  calculated  for  the  case  of 
normal-distributed  loads  as 


Figure  2:  Comparison  of  the  number  of  cycles  in 
DLDR  phases  I  and  II  to  LDR  cycles  to  failure  - 
total  DLDR  phase  I  and  II  cycles  at  any  given  load 
level  equals  LDR  cycles  but  high  values  of  S/E  in 
phase  I  limit  lower  values  in  phase  II. 


Figure  2  further  illustrates  the  DLDR  concept  where 
high  values  of  S/E  during  phase  I  would  limit  the 
number  of  cycles  to  failure  for  lower  values  of  S/E 
during  phase  II.  However,  for  simple  cases  involving 
a  single  load  level,  the  number  of  cycles  is  split 
between  the  two  phases  and  the  total  number  of 
cycles  remains  unchanged. 

DOUBLE-LINEAR  RELIABILITY  THEORY 

This  paper  aims  to  establish  a  double-linear  reliability 
theory  based  on  expressions  using  generalized 
a  -functions  or  similar  for  each  DLDR  phase.  The 
inverse  of  phase  I  cycles  for  the  t-th  cycle  applied  to 
the  y'-th  component  may  be  expressed  as 


where  m0  =  —mQ.  In  accordance  with  Manson- 
Halford  DLDR  theory,  the  y'-th  component  is 
predicted  to  enter  phase  II  after  accumulating  phase  I 
damage  of  D,j  =  1. 


The  inverse  of  phase  II  cycles  for  the  i-th  cycle 
applied  to  the  y- th  component  may  be  expressed  as 


^  0  otherwise 

for  which  a  new  class  of  generalized  DLDR-based 
phase  II  a-functions  is  defined  herein  for  v  <  z0  as 


a*DLDRn(z0,m,v,y,m0)  = 

(z  —  vy 


r-  W. 

■'zo 


1  —  exp[y(z  —  u)m°] 


0(z)dz 


(19) 


^  0  otherwise 

for  which  generalized  a-functions  are  not  applicable 
without  first  incorporating  the  effects  of  the 
denominator  terms  absent  from  similar  expressions 


As  in  the  case  of  phase  I,  generalized  DLDR-based 
phase  II  a-functions  may  be  applied  to  phase  II 
damage  calculations,  as  follows 

Dn,j  —  n n  9e,ii  (^j)  (20) 


where  Dn  j  is  the  Manson-Halford  DLDR  phase  II 
damage  fraction  and  the  corresponding  phase  II 
damage  rate  (?£•,//(£} )  is  calculated  for  the  case  of 
normal-distributed  loads  as 


Per  Manson-Halford  DLDR  theory,  the  y'-th 
component  is  predicted  to  fail  after  entering  phase  II 
and  accumulating  phase  II  damage  of  DUj  —  1. 


It  is  recommended  that  the  integral  for  %  be 
evaluated  using  a  Newton-Coats  formula,  such  as  a 
1000-step  4th-order  formula  ( i. e. .  Boole’s  method,  see 
[13]  for  details). 


Implementation  of  Gauss-Laguerre  requires  a  change 
of  coordinates.  As  was  the  case  in  [7],  let  x  =  |z2 
and  xshift  —  rzx'1  for  the  generalized  DLDR 
a-function.  One  may  evaluate  a2  as 


p  00 

a2  =  f{x ) exp(— 


x)  dx 


xshift 


(26) 


Zn 

Wi  f(xt  +  xshift) 

1=1 


For  completeness,  similar  DLDR-based  functions  for 
use  with  Weibull  loads  (extending  the  b-fiinction 
tools  from  [5]  and  [7])  are  developed  in  Appendix  A. 

EVALUATION  TECHNIQUES 


where 


fix)  = 


1 

yjAnx 


(V2x  —  v)m 


(27) 


Although  the  author  introduced  analytic  evaluation  of 
a  certain  subset  of  a-functions  in  [5],  the  author 
recommends  numerical  integration  techniques 
introduced  in  [7]  be  applied  for  evaluation  of  the 
generalized  DLDR  a-functions,  a*DLD mi*, •,»,•,•)  and 

aDLDRll  (•<•<•<•<•)■ 


and  values  for  wt  and  xt  are  provided  in  [14]. 

Evaluation  of  generalized  DLDR  phase  II 
a-functions,  cl*dldrii  (•,«,•,*,•)  may  be  accomplished 
with  similar  methods  after  making  the  following 
modifications  to  certain  of  the  above  expressions: 


In  [7],  evaluation  of  a-functions  with  non-integer  m 
values  is  accomplished  using  15-point  Gauss- 
Laguerre  quadrature  formulas  (see  [13]  and  [14]  for 
details).  To  correct  for  inaccuracies  in  use  of  Gauss- 
Laguerre  for  evaluating  the  a-function  for  small 
values  of  \z0\  and  m,  [7]  recommends  that  the 
a-function  integral  be  evaluated  in  two  steps. 
Adapting  the  methods  of  [7]  to  the  generalized 
DLDR  phase  I  a-functions  integral  results  in  the 
following: 

a*DLDRi  0?o.  m,  v,  y,  m0)  =  ax  +  a2  (22) 

where 

f  zi  ( z  —  v)m 

ai  JZo  explyiz  -  v)m°]  (23) 

az  =  i  exp[y(z  —  v)m°]  ^z)dz  W 
and 

Zi  =  1  +  +  lm  +  \z02  (25) 


a*DLDRii  (z0-  m,  v,  y,m0)  =  at  +  a2  (28) 
rzi  ( z  -v)m 

ai  —  I  - ^ — r  4>iz)dz  (29) 

JZq  1  —  exp[y(z  —  v)m°]  v 

[°°  (z-v)m 

=  l  1  -  exp[y(z  -  P)"°]  <301 


fix) 


1  (V2x  —  v)m 

V47TX  i  —  exp  [y(V2x  —  v )  °j 


(31) 


Numerical  procedures  and  all  other  expressions 
remain  unchanged  from  those  used  in  evaluation  of 
the  generalized  DLDR  phase  I  a-function. 


APPLICATION  NOTES 


In  this  section,  generalized  DLDR-based  phase  I  and 
phase  II  a-functions  are  applied  to  three  special  cases 
of  fatigue  reliability  problems.  Each  case  adapts 
results  in  [7]  to  study  double-linear  reliability  in  the 
presence  of  load  and  strength  variation.  Two  cases  of 
deterministic  usage  are  considered,  one  for  usage 
characterized  by  a  single  load  condition  and  the 
second  for  a  distributed  usage  over  multiple  load 


conditions.  Also,  a  probabilistic  usage  case  is 
considered  for  a  single  load  condition. 

An  approach  to  apply  LDR-based  a-function  tools  to 
multivariate  probabilistic  usage  problems  is  presented 
in  [7].  However,  due  to  the  complexity  of  the 
effective  composite  DLDR  damage  rate  for  multiple 
load  conditions,  a  similar  approach  is  not  available 
for  application  to  the  DLDR  case  at  this  time. 


Single  Load  Condition  Deterministic  Usage 

Application  is  first  considered  for  deterministic  usage 
with  a  single  load  condition.  At  component  failure, 
accumulated  damage  D(  in  the  j-th  component  is 
equal  to  some  threshold  value,  S  (typically  failure  is 
considered  to  occur  at  S  —  1).  At  failure,  the  phase  I 
and  phase  II  damage  fractions  are  also  at  the 
threshold,  namely:  Z);y  =  S  and  Du j  —  S. 


The  reader  may  object  to  the  concept  of  using  the 
same  threshold  value  for  phase  I  and  II.  In  fact,  the 
sequence  of  damage  accumulation  is  that  phase  I  is 
first  accomplished  prior  to  initiation  and  completion 
of  phase  II.  However,  the  author’s  concept  of  the 
threshold  is  such  that  both  phases  have  been 
completed  at  failure  by  reaching  the  threshold  (each 
in  turn).  Sequence  is  only  important  if  the 
deterministic  usage  changes  with  time  and  varies 
between  phase  I  and  phase  II.  In  such  cases,  the 
phase  I  and  II  damage  rates  should  be  calculated 
according  to  the  appropriate  usage  spectrum  during 
each  phase  and  the  following  discussion  would 
remain  applicable. 


The  number  of  cycles  to  failure  (specifically,  to 
achieve  Dj  =  S )  is  calculated  as 


n  =  nI  +  nu 


(32) 


Let  E  be  randomly  distributed  with  density  function 
fE(E).  The  probability  that  the  damage  in  a  given 
component  is  less  than  the  failure  threshold  S  is 

Pr(D<5)  =  FD(5)=  j  fE(E)dE  (35) 

All  values 
of  E  where 
gEn{E,n)<S 


where  Pr(£>  <  S )  is  the  component’s  reliability,  R. 


Conversely,  the  probability  that  the  damage  in  a 
given  component  exceeds  the  failure  threshold  S  is 

Pr(D  >  S')  —  j  fE(E)dE  (36) 

All  values 
of  E  where 
gEn{E,n)>8 


where  Pr(£>  >  S )  is  the  probability  of  component 
failure.  Considering  a  specific  critical  endurance 
limit  Es  such  that  gEn(Es,ri)  =  S  for  a  given  value 
of  n,  one  may  clarify  the  integration  as 


Es 

Pr(D  >  S)  —  j  fE(E)dE  (37) 

—  CO 

which  is 

Pr (D  >  S)  =  FE (Eg)  (38) 

by  definition,  where  FE(»)  is  the  cumulative 
endurance  limit  distribution.  It  is  important  to  clarify 
that  the  particular  type  of  fatigue  strength  distribution 
is  not  specified  (i.e.,  the  methods  of  this  paper  are 
applicable  to  normal  strength,  log-normal  strength, 
Weibull  strength,  etc.). 

Single  Load  Condition  Probabilistic  Usage 


Solving  for  S  yields  the  accumulated  damage  as 

Dj  =  gEn(Dj,n )  =  ngE(Ej )  (33) 


where 


9e(Ej) 


9E.i{Ej)9  E,II  C Ej ) 
9E,i{Ej)  +  gE,u{Ej) 


(34) 


With  probabilistic  usage,  the  number  of  cycles  of 
each  load  condition  at  retirement  is  also  a  random 
variable.  The  simple  case  of  a  single  damaging  load 
condition  is  considered,  where  the  accumulated 
damage  may  be  expressed  as 

Dj  =  njgE  (Ej)  =  gEn  ( E} ,  ny)  (39) 


9e  (Ej) 


9E,i(Ej)gE  ,n(Ej) 
9E,i(Ej)  +  gE,u(Ej) 


may  be  considered  the  effective  composite  DLDR 
damage  rate. 


where 


(40) 


is  the  effective  composite  DLDR  damage  rate. 

For  this  case,  one  may  substitute  the  effective 
composite  DLDR  damage  rate  directly  into  the 
discussion  of  probabilistic  usage  in  [7],  where  the 
probability  that  the  damage  is  above  some  threshold 
S  (probability  of  component  failure)  is  calculated  as 


is  the  effective  composite  DLDR  damage  rate  for 
multiple  load  conditions. 

As  before,  one  may  consider  a  specific  critical 
endurance  limit  Es  such  that  gEn(Es,n)  =  S  for  a 
given  value  of  n  and  solve  for  the  probability  of 
component  failure  as 


Pr(D  >  5)  =  Fb(0) 

n  CO 

+  /*(£) 
^0 


Pr  {D>8-)  =  Fe{Es)  (48) 

where  FE(»)  is  the  cumulative  strength  distribution. 


where  E  and  n  are  independent  and  randomly 
distributed  with  density  functions  fE(E)  and  fn(n) 
and  cumulative  distributions  Fe{E )  and  Fn(n), 
respectively. 

Multiple  Load  Conditions  Deterministic  Usage 

For  the  case  of  deterministic  usage  with  p  multiple 
(potentially  damaging)  load  conditions  in  the  usage 
spectrum,  the  Manson-Flalford  DLDR  phase  I  and 
phase  II  damage  fractions  are  calculated  as  follows: 


EXAMPLES 

Examples  are  used  demonstrate  application  to  single 
and  multiple  load  cases. 

Tables  1  through  4  present  an  example  intended  to 
allow  demonstration  of  the  double  linear  reliability 
methods  developed  in  this  paper.  Table  1  presents 
the  problem  parameters  where  k  and  m  correspond  to 
the  curve  shape  shown  in  Figure  2  and  the  mean  load 
is  twice  the  mean  endurance  limit. 


D,,j  =  n,  lT  gJJ(Ej)  (42) 


and 

Dn,j  —  tin  9e,ii{E]) 


respectively,  where 


f  = 


rfii 

f2 


LfpJ 


(43) 


(44) 


is  a  vector  of  load  condition  weighting  factors 
describing  the  usage  spectrum,  and  the  vectors 
gE'i{Ej )  and  gE  l}{Ej)  contain  corresponding  DLDR 
phase  I  and  phase  II  damage  rates  for  each  load 
condition  in  the  usage  spectrum. 


The  number  of  cycles  to  failure  (Dj  —  S)  is 


(45) 


n  =  n,  +  nn  =  S 


■  +  ■ 


KirgiftEj)}  ~gEJi{Ej)}) 
Solving  for  S  yields  the  accumulated  damage  as 

Dj  =  gEn{Ej,n)  =  ngE(Ej )  (46) 


where 


9e(Ej)  = 


[iT  9E,i(Dj)}  9e, //(£))} 

r?T - .  (?t - 


For  those  unfamiliar  with  a-functions  presented  in 
[5]  and  [7],  a  baseline  solution  is  developed  via 
application  of  traditional  methods  with  results 
provided  in  Table  2.  A  fatigue  life  of  just  over 
90,000  cycles  is  calculated  using  a  fatigue  curve  with 
mean  minus  three  sigma  working  endurance  limit, 

Ew,  and  a  load  level  corresponding  to  95th  percentile 
(or  top  of  scatter)  vibratory  loads,  Sg5. 

Table  3  presents  results  for  the  example  problem 
using  double  linear  theories  developed  in  this  paper. 

A  fatigue  life  of  nearly  68,000  cycles  corresponds  to 
“six  nines”  of  reliability  ( R  =  0.999999)  for  the 
problem  characterized  by  parameters  in  Table  1, 
which  is  25%  less  than  the  fatigue  life  predicted 
using  traditional  methods  (see  Table  2).  The  reader 
should  attempt  replicating  the  single  load  case  results 
provided  in  Table  3  prior  to  attempting  more  realistic 
problems  involving  multiple-load  cases. 

Table  4  presents  results  for  LDR  methods  developed 
in  [7].  These  results  are  presented  to  allow  the  reader 
to  test  the  ability  to  duplicate  simpler  LDR  results 
prior  to  attempting  DLDR  methods.  As  one  might 
expect,  comparison  of  DLDR  and  LDR  results  in 
Tables  3  and  4,  respectively,  indicates  little 
difference  between  DLDR-based  and  LDR-based 
reliability  methods  when  applied  to  problems  with  a 
single  damaging  load  condition. 


(47) 


Table  1:  Parameters  Used  for  a 
Single  Load  Condition  Example 


Table  3:  Single  Load  Condition  Example  Results 
with  a  Target  Reliability  of  “Six  Nines”  (0.999999) 


Parameter 

Value 

Variable  Name 

Value 

k 

1.0  X  106 

D1-1, 

0.0622398 

m 

2 

P 

-62.7717 

P-E 

5 

Q 

-0.451408 

COVE 

0.1 

Ecrit 

2.62329 

Ms 

10 

z0 

-7.37671 

COVs 

0.1 

V 

-10.0000 

^ref 

1.0  x  103 

y 

-0.0514247 

N* 

1.0  x  106 

m0 

0.902817 

a*DLDRi  (zo>  m,v,y,  m0) 

153.562 

Table  2:  Results  for  Traditional  LDR  Method' 

a*DLDRii(zo,m,  v,y,m0 ) 

297.367 

Variable  Name 

Value* 

9e,i 

2.23147  X  10-5 

Ew 

3.50000 

9e,ii 

4.32117  X  10-5 

•S95 

11.6449 

9e 

1.47155  X  10-5 

n 

90337 

n 

67955 

(per  traditional  method) 

1  Calculation  is  based  on  mean  minus  3  sigma  working 
fatigue  strength  curve  and  95-th  percentile  top  of  scatter 
vibratory  load. 

*  Note:  intermediate  results  in  this  and  each  of  the  other 

Table  4:  Comparative  Results  for  LDR  Case  with 
Six  Nines  Target  Reliability 

tables  in  this  paper  are  rounded  to  6  digits  and  calculated 
fatigue  lives  are  rounded  to  the  nearest  whole  cycle. 

Variable  Name* 

Value 

a*(z.0,m,  v ) 

101.000 

9e 

1.46767  X  10-5 

n 

68135 

*  Note,  variables  in  this  table  are  defined  per  methods 
presented  in  [7]  rather  than  double  linear  equations  herein. 


To  demonstrate  application  to  problems  involving 
multiple  load  conditions,  an  example  is  borrowed 
from  [9]  where  four  load  conditions  are  applied 
which  would  produce  nominal  lives  of  103,  104,  10s, 
and  106  cycles  (for  mean  loads  and  strength).  As  in 
[9],  each  block  of  loads  contains  1  percent  of  the 
nominal  failure  cycles  for  each  loading  level 
individually.  Tables  5  and  6  adapt  this  example  to 
demonstrate  the  methods  developed  in  this  paper.  It 
is  noted  that  the  author  has  modified  the  problem  by 
assuming  that  each  “block”  of  loads  corresponds  to 
400  flight  hours  (selected  to  correspond  to  100  hours 
per  damaging  load  condition). 

Results  are  presented  in  Figure  3,  where  four 
solutions  are  shown  for  each  of  the  LDR  and  the 
DLDR  theories: 

•  Nominal  method  (mean  endurance  limit) 

•  Traditional  method  (working  endurance  limit) 

•  Probabilistic  method  (a-functions) 

•  Probabilistic  method  (a-functions)  with 
additional  load  below  working  endurance  limit 

First,  a  life  is  calculated  using  nominal  methods 
applied  for  the  mean  endurance  limit  (100  load  units), 
resulting  in  a  10,000  hour  life  for  LDR  and  a 
4,596  hour  life  for  DLDR.  The  LDR  result  is  as 
expected  when  considering  that  four  load  conditions 
are  applied  per  400  hour  block  which  each  produce 
1%  damage  fraction  per  block  (4%  damage  per  400 
hours  ratios  to  100%  damage  per  10,000  hours). 

Also,  the  DLDR  result  corresponds  to  within  0.2%  of 
that  presented  in  [9]  (after  adjusting  for  400  hours  per 
load  block). 

Second,  using  the  working  endurance  limit  (70  load 
units),  traditional  methods  result  in  lives  of  4,900 
hours  and  2,089  hours  for  LDR  and  DLDR  cases, 
respectively. 

Next,  the  methods  provided  in  [7]  and  herein  are 
applied  to  the  problem  based  on  the  strength 
Coefficient  of  Variation  (COV)  and  the  load  COV 
specified  in  Table  5.  To  maintain  six  nines  of 
reliability,  the  resulting  lives  are  2,725  hours  and 
1,222  hours  for  LDR  and  DLDR  cases,  respectively. 

Finally,  the  usage  spectrum  in  Table  7  is  considered, 
where  an  additional  load  is  included  with  load  level 
below  the  working  endurance  limit,  Ew,  and 
occurring  10,000  cycles  per  400  hour  block.  This 
level  of  load  would  be  neglected  by  traditional 
methodologies  where  loads  below  Ew  do  not  impact 
the  calculation.  However,  the  probabilistic  methods 
provided  in  [7]  and  herein  are  impacted  at  high 


Table  5:  Parameters  Used  for  a 

Multiple  Load  Condition  Example 

Parameter 

Value 

k 

1.0  X  106 

m 

2 

Me 

100 

COVE 

0.1 

Es 

varies 

COVs 

0.1 

^ref 

1.0  x  103 

N* 

1.0  x  106 

Table  6:  Usage  Spectrum  for  Example  with 
Multiple  Load  Conditions 

Index,  p 

Mean  Load 
Level 

Nominal 
Cycles  to 
Failure 

Cycles  per 
400  hour 
Block 

1 

3162.28 

103 

10 

2 

1000.00 

104 

100 

3 

316.228 

10s 

1000 

4 

100.000 

106 

10000 

Table  7:  Alternative  Usage  Spectrum  for 
Example  with  Multiple  Load  Conditions 


Index,  p 

Mean  Load 
Level 

Nominal 
Cycles  to 
Failure 

Cycles  per 
400  hour 
Block 

1 

3162.28 

103 

10 

2 

1000.00 

104 

100 

3 

316.228 

10s 

1000 

4 

100.000 

106 

10000 

5 

66.6667 

00 

10000 

DLDR-based  reliability 
LDR-based  reliability 


DLDR,  with  additional  load  <  (p-3o)E 
LDR,  with  additional  load  <  (p-3a)E 


Life  (hours) 


Figure  3:  Results  for  the  Multiple-Load-Condition  Example  (Adapted  from  [9])  with  400  Hour  Load  Blocks 


reliability  levels  based  on  overlap  of  the  load  and 
strength  distribution  tails.  To  maintain  six  nines  of 
reliability,  the  resulting  lives  are  2,455  hours  and  934 
hours  for  LDR  and  DLDR  cases,  respectively.  The 
impact  of  the  additional  load  level  below  Ew  is 
magnified  in  the  DLDR  case  due  to  reductions  in 
cycles  to  failure  at  low  load-levels  during  phase  II 
(see  Figure  2). 

SENSITIVITY  STUDIES 

This  section  illustrates  the  sensitivity  of  the  DLDR 
cumulative -damage  reliability  tools  to  load  and 
strength  variations.  To  avoid  complex  sensitivity 
studies  based  on  multiple  load  condition  interactions, 
the  single  load  condition  example  with  parameter 
values  provided  in  Table  1  is  taken  as  a  baseline  for 
comparison.  Although  single  load  condition 
examples  do  not  take  full  advantage  of  DLDR 
theories,  these  sensitivity  studies  fully  apply  the 
DLDR  method  and  may  be  used  to  understand  the 
impact  of  changes  in  any  given  load  condition. 


Figures  4  and  5  demonstrate  the  increase  in 
probability  of  failure  (decrease  in  reliability)  as  the 
number  of  cycles  is  increased  and  with  increasing 
load  ratio.  Although  the  nominal  ASTM  log-log 
strength  curve  demonstrates  a  sharp  comer  at 
Nel  —  10c,  the  author  notes  that  probabilistic 
methods  tend  to  smooth  the  corner. 

Figure  6  presents  results  for  increases  in  load  COV 
with  six  nines  reliability.  Load  variability  effects  are 
more  significant  below  the  endurance  limit  where  the 
curves  in  Figure  6  “fan  out”  with  additional  cycles 
and  increased  load  variation.  Figure  6  demonstrates 
the  importance  of  understanding  load  distributions 
for  sustained  steady-state  conditions  such  as  level 
flight,  which  are  typically  underrepresented  in  a 
composite  worst  case  usage  spectrum.  Figure  7 
presents  results  for  increases  in  strength  COV. 
Comparison  of  the  14%  strength  COV  case  from 
Figure  7  to  the  50%  load  COV  case  from  Figure  6 
indicates  greater  sensitivity  to  strength  variation  at 
each  load  level  (see  [5]  for  similar  LDR  results). 


HS/fiE  =  Probability  of  Failure 


yiS/jaE  =  8  - fiS/ftE  =  4  . yiS/yiE  =  2 

- yiS/jaE  =1  - jiS/jiE  =  0.67 - yiS/yiE  =  0.5 


n 

4:  Probability  of  Failure  as  a  Function  of  Cycles  for  Various  Ratios  of  Mean  Load  to  Mean  Strength, 
HS/ nE,  with  k  =  106,  m  =  2,  COVE  =  0. 1,  COVs  =  0. 1,  Vre/  =  103,  and  N*  =  106 


Mean  Fatigue  Curve - 4  nines  - 5  nines 

6  nines  - 7  nines  - 8  nines 


n 


Figure  5:  Ratio  of  Mean  Load  to  Mean  Strength,  fis/fiE,  as  a  Function  of  Cycles  for  Various  Reliability 
Levels,  with  k  =  106,  m  =  2,  CW£  =  0. 1,  C0Fx  =  0. 1,  Wre/  =  103,  and  N *  =  106 


3ri/sr1  =  Hri/sr/ 


Mean  Fatigue  Curve - CO  VS  =  0.1 

COYS  =  0.3  - COYS  =  0.4 


- COVS  =  0.2 

- COVS  =  0.5 


n 


6:  Ratio  of  Mean  Load  to  Mean  Strength,  j us/fiE,  as  a  Function  of  Cycles  for  Various  Levels  of  Load 
Variation,  with  six  nines  reliability,  k  —  106,  m  —  2,  COVE  —  0. 1,  Nref  —  103,  and  JV*  =  106 


Mean  Fatigue  Curve - COVE=0. 1  - COVE=0. 1 1 

COVE=0.12  - COVE=0.13  - COVE=0.14 


n 


Figure  7:  Ratio  of  Mean  Load  to  Mean  Strength,  fis/ nE,  as  a  Function  of  Cycles  for  Various  Levels  Strength 
Variation,  with  six  nines  reliability,  k  —  106,  m  —  2,  COVs  —  0. 1,  Nref  —  103,  and  N*  —  106 


CONCLUSIONS 


The  author’s  conclusions  are  as  follows: 

1 .  The  work  described  in  this  paper  extends  the 
generalized  a-function  tools  developed  in  [5]  and 

[7]  to  allow  calculation  of  fatigue  reliability  based 
on  the  Manson-Halford  DLDR,  thus  formulating  a 
new  double-linear  reliability  theory. 

2.  A  sensitivity  study  has  been  performed  which 
indicates  the  importance  of  understanding  the 
strength  distribution,  especially  the  strength  COV. 

3.  A  multiple-load-condition  example  clearly 
demonstrates  that  the  selection  of  cumulative 
damage  theory  impacts  reliability. 

a.  Probabilistic  DLDR  methods  demonstrate 
significant  reductions  in  fatigue  life  to 
maintain  six  nines  of  reliability  when 
compared  to  LDR  methods. 

b.  The  importance  of  understanding  and  applying 
the  appropriate  LDR,  DLDR,  or  other  material 
characterization  is  apparent. 

4.  Future  efforts  should  include  the  following: 

a.  The  potential  benefits  of  re-introducing  £>/_// 
spectrum  dependency  into  more  complex 
methods  should  be  assessed. 

b.  The  applicability  of  double-linear  reliability 
methods  described  in  this  paper  would  be 
significantly  improved  by  the  development  of 
an  approach  to  apply  DLDR-based  methods  to 
multivariate  probabilistic  usage  problems. 
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APPENDIX  A 


with 


This  appendix  presents  double-linear  reliability 
theory  for  use  with  Weibull  distributed  loads  by 
extending  the  6-function  methods  of  [5]  and  [7]. 
Although  the  material  in  this  appendix  is  not  required 
to  understand  and  apply  the  a-function  methods  to 
problems  involving  normally  distributed  loads,  these 
tools  are  provided  for  completeness. 

Additional  Notation  used  in  Appendix  A: 

6  (•,«,•)  Benton  6-function  defined  in  [5] 

6*(v •,»,•)  generalized  6-function  defined  in  [7] 
^dldri (•/•/•/•/•>•)  generalized  DLDR-based  phase  I 
6-function  defined  herein 

b*DLDMi  (•/•/•/•.*/•)  generalized  DLDR-based  phase  II 
b -function  defined  herein 
u  generalized  6-function  parameter 

w  “standardized”  Weibull  load  amplitude 

W{  weighting  value  used  in  Gauss-Laguerre 

quadrature 

w0  “standardized”  Weibull  damaging  load 

amplitude 

P  Weibull  slope  (for  load  amplitudes) 

y(»,«)  “standardized”  Weibull  probability 

density  function 

S  Weibull  minimum  expected  value  (for 

load  amplitudes)  -  also  used  above  as 
accumulated  fatigue  damage 
r]  Weibull  characteristic  value  (for  load 

amplitudes) 


(A-3) 

ip  wP  1  exp(—  w^)  if  w  >  0 
y(w, /?)  =  j 

vO  otherwise 


Similarly,  the  DLDR phase  II  damage  rate,  gEn(Ej), 
is  calculated  for  the  case  of  Weibull  loads  as 


where  the  new  class  of  generalized  DLDR-based 
phase  II  6-functions  is  defined  for  u  <  max(w0, 0)  as 

Kldrii  Oo,  m>  P>  u,  y,  m0)  =  (A-5) 

r  O  -  u)m  ,  N 

I  - ^ — r  Y(w,P)dw 

Jmax(w0,0)  1  -  exp[y(w  -  U)m o] 

It  is  noted  that  the  DLDR  phase  I  and  phase  II 
damage  rates  for  Weibull  loads,  gE[{Ej)  and 
gE,n{Ej),  respectively,  may  be  substituted  into 
equations  used  for  various  special  cases  in  the 
Application  Notes  section  of  the  main  body  of  this 
paper  to  calculate  the  effective  composite  DLDR 
damage  rate,  gE{Ej). 


Double-Linear  Reliability  Theory  for  Weibull  Loads 

The  DLDR  phase  I  damage  rate,  gE  I  (£)■),  is 
calculated  for  the  case  of  Weibull  loads  as 


where  and  P,k,  Q,  and  m0  are  as  defined  in  the  main 
body  of  this  paper;  /?,  g,  and,  S  are  Weibull 
distribution  parameters;  and  the  new  class  of 
generalized  DLDR-based  phase  I  6 -functions  is 
defined  for  u  <  max(w0, 0)  as 

Kldri  Oo,  m,  p,u,y,  m0)  =  (A"2' 

r°°  (w  -  u)m 

- ^ — ^y{w,p)dw 

-'max(w0,0)  6Xp[y(w  —  u)m°] 


